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Asymptotic  Properties  of  a  Sensor  Ailocation  Modei 


D.  P.  Gaver 
P.  A.  Jacobs 
M.  Youngren 

0.  Introduction  and  Summary 

Sensor  resources  are  to  be  allocated  among  /  non-overlapping  but  possibly 
contiguous  geographical  locations  called  nodes.  Each  sensor  look  at  a  node  results 
in  an  observation,  with  error,  of  the  number  of  units  that  are  at  the  node.  In 
military  applications  units  may  be  individual  assets  such  as  tanks  or  ships,  or 
possibly  small  groups  organized  as  platoons  or  companies  and  in  geographical 
proximity;  in  ecological  applications,  they  might  be  animals,  singly  or  in  groups, 
or  particular  vegetation  types.  The  observations  are  then  to  be  used  to  estimate 
the  numbers  of  units  at  each  node.  The  problem  is  to  allocate  sensor  resources  so 
as  to  minimize  a  measure  of  the  error  of  the  estimates,  and  of  the  estimate  of  the 
sum  of  the  numbers  at  all  nodes.  It  is  assumed  that  the  units  remain  on  the  nodes 
and  do  not  migrate,  although  for  some  sensor  types  the  units  must  be  in  motion 
on  the  node  at  least  some  of  the  time  in  order  for  detection  to  take  place. 

The  objective  of  this  paper  is  to  propose  and  study  adaptive  allocation  of  sensor 
effort  in  such  a  way  as  to  focus  the  sensor's  attention  sequentially  and 
purposefully  on  nodes  so  as  to  pay  most  attention  to  those  nodes  about  which 
die  greatest  uncertainty,  or  interest,  currently  prevails.  This  problem  has  features 
in  common  with  adaptive  bin-packing  problems  (cf.  Gaver  et  ah,  1995)  and 
adaptive  allocation  of  customers  to  servers  (cf.  Gaver  et  al.,  1993). 


Numerical  examples  show  that  a  properly  selected  sequentially  adaptive  rule 
will  provide  estimates  of  improved  precision. 

A  related  investigation  that  has  recently  come  to  our  attention  is  that  by 
Thompson  and  Seber  (1994).  Their  procedures  can  also  be  analyzed  using  the 
methodology  proposed  in  this  paper. 

1.  The  Model 

There  are  I  nodes.  Assume  that  node  i  contains  r,-  units.  Sensor  resources  are 
allocated  to  one  of  the  I  nodes  at  times  that  occur  according  to  a  Poisson  process 
with  rate  A.  The  sensor  resource  is  allocated  to  node  i  with  probability  af  which  is 
tailored  to  depend  on  past  allocations  in  a  purposeful  way.  Let  Z„(i)  be  the 
observation  of  node  i;  a  simple  model  is  that  {Zn(i);  n  =  1, 2, ...}  are  iid  with 
binomial  distribution  with  ri  trials  and  known  probability  of  success  pi  which  is 
the  probability  of  unit  detection  at  node  i.  We  assume  for  illustration  that  the 
number  of  imits  on  a  node  does  not  change,  although  an  adaptive  scheme  of  the 
type  proposed  should  effectively  follow  changes  in  node  population. 

Let  Ni(f)  be  the  number  of  times  in  [0,/]  that  node  i  is  observed  by  a  sensor. 
Let 

V(i;t)  =  Zi(i)+ ...  + 

be  the  sum  of  all  the  observations  during  [0,f]  for  node  i. 

Under  some  conditions  Vii;  t)  has  a  binomial  distribution  with  Niit)ri  trials 
and  probability  of  success  pi,  given  Niit),  although  this  assumption  may  not  be 
especially  accurate  in  general.  Assuming  it  to  be  adequate  for  the  moment,  an 
estimate  for  the  number  of  units  at  node  i  based  on  observations  made  during 


with 


=  (1.2) 

Niit)  pf 

_  _  fiit){i-pi) _  v{ut){i-pi) 

~  Niitfpf 

There  has  been  much  statistical  attention  paid  to  estimating  the  number  of  trials 
in  a  binomial  distribution,  given  probability  of  success;  for  recent  discussion  see 
Hall  (1994).  Here  we  use  the  simplest  such  estimator. 

Now  consider  the  following  adaptive  allocation  rule.  A  sensor  resource  arriving 
at  time  t  is  allocated  to  node  i  with  probability 

(1.4) 

y=i 

where  the  hi  are  strictly  positive  sufficiently  smooth  functions  such  that 
hjicPx,  cVy)  =  cPhjix,  y). 

One  possible  form  for  hj  is 


hi{x,v)=  — +  for 

X  Pi  . 


y>0  and  ai>0; 


if  it  is  assumed  that  given  N,(0,  Vjit)  is  binomially  distributed  with  mean 
then  from  (1.3)  it  follows  that  for  this  function  hi,  sensor  resources  tend  to  be 
allocated  to  those  nodes  for  which  the  variance,  and  hence  its  square  root,  the 
standard  error  of  the  estimated  number  of  units,  is  the  largest;  the  probability  of 
allocation  to  the  most  uncertain  node  (by  this  measure)  increases  rapidly, 
approaching  unity  as  yinaeases.  This  tends  to  bring  down  that  standard  error 
quickly,  and  to  equate  standard  errors  of  the  estimates  across  nodes.  Clearly, 
alternative  measures  of  overall  sensor  performance  are  feasible,  and  possibly 
desirable,  such  as  ones  that  endeavor  to  equalize  fractional  or  percent  error  on 
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nodes,  or  ones  that  also  respond  to  an  independent  measure  of  importance  of  the 
units  on  a  node.  Additionally,  node  contents  may  be  of  various  types,  which  can 
be  considered.  Thus,  the  present  discussion  is  of  an  illustration  of  an  adaptive 
allocation  scheme. 

The  purposeful  allocation  of  (1.4)  introduces  dependence  between  {N,(t);  t  >  0) 
and  {Vi(t);  t  ^  0}.  Asymptotic  results  for  the  means  E[Ntit)]  and  E[Vi(0]  as  the  rate 
of  the  Poisson  process  A  ^  are  obtained  in  Section  2.  It  is  shown  that  the 
purposeful  allocation  estimate  fi{t)  is  asymptotically  imbiased.  Section  3  presents 
results  for  the  asymptotic  means  for  specific  form  of  function  hi  (1.5).  Section  4 
discusses  asymptotic  results  for  the  second  moments  of  {N,(f);  f  ^  0}  and  {VKO; 

im  ■ 

i 

Section  5  describes  a  simpler  Poisson  approximation.  Section  6  presents  results 
from  simulation  experiments. 


f  ^  0}  and  presents  approximate  expressions  for  the  and  Var 


2.  First-Moment  Calculations  and  Asymptotics 

Note  that  for  >  0 

£[N(((  +  ft)|N{(),  V(f)]  =  N,(f)  +  XoiiN,  V)h . 

Thus 

E[N,(f  +  h)]  =  E[N,(f)]  +  XE[ai{N,V)h]. 
Assuming  derivatives  exist,  we  have 

l£[Nj(()]  =  A£[oi(lV(/),V(())]. 

Similarly, 

E[v;(f + hp{tlv{t)] = v;(f) + ^{N,V)riPi. 


(2.1) 


(2.2) 


(2.3) 

(2.4) 
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Thus, 


E[v;-(f  +  h)]  =  E[v;(0]  +  XhriPiE[ai{N,V)h].  (2.5) 

Assuming  derivatives  exist,  we  have 

|E[V((()]  =  riPi^h(N,F)].  (2.6) 

Assume  oq  is  of  the  form  (1.4)  where  hi  is  sufficiently  smooth  with 
hiicPx,  cVy)  =  cVhiix,  y)  for  constants  c  and  hj^x,  y)  >  0  for  all  x,  y  >  0.  We  have 


Ift/KW-VyW) 

/ 

ft,(Ni(()/A,v;(0A) 


Assume  lim  Nj(t]lX  =  lim  Vi{t)/X  =  »,■{().  lini 

X—>oo  X-^oodt  at 

and  ^Im  ^E[v;(f)/A]  =  ^ i;,(f) . 

Dividing  both  sides  of  (2.3)  and  (2.6)  by  A  and  letting  A  the  bounded 
convergence  theorem  yields 

(2.8, 


±v.(t)=rp . ■ 

dt  ‘  ‘  ^  ^hj{mj{t),Vj{t)) 


Thus, 


Vi{t)  =  riPimi{ty 


(2.10) 
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3.  Purposeful  Allocation 

The  equations  for  (2.8),  can  be  solved  for  special  functions  hi.  For 


example,  assume 


Since 


hi{x,y)=  — +  o  7  ^  fory^O  andfl,  >0. 
^  Pi  . 


Vi{t)  =  mi{t)Tipi, 

■'  ''  mj(()  I 


In  this  case  (2.8)  can  be  rewritten  as 


^^•(0  _  ai  +  rjCj  ^ji^) 
^my(f)  [aj  +  rjCjmiit) 


where  c,  =  (l-p,)/p,-. 
The  functions 


ii{t)  =  K{ai  +  riCiY^^'^^^h 


with 


are  a  solution  to  (2.8). 

Thus,  in  the  limit  as  f  -» «>,  the  probability  that  a  sensor  resource  at  time  t 


looks  at  node  i  is 


V/  ^  \r/(r+i)’ 
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If  7=  0,  then 


«,•  =  !// 


and  the  allocation  is  equally  likely. 
If  7-4  00,  then 


If,  further,  the  probability  of  detection  on  node  i,  pi,  is  constant  for  all  the  nodes, 
then  the  probabilistic  allocation  for  7-^ » is  roughly  proportional  to  the  number 
of  units  on  the  node. 

4.  Scaling  and  Approximate  Variances  of  Estimators 


Xy(() - (4.1) 

1^'  (4.2) 

In  Appendix  A,  moment  generating  functions  are  used  to  show  the 
asymptotic  normality  of  {(Xy(f),  V;(0);  f  ^  0}  as  A  -4  00.  Further,  differential 
equations  for  their  second  moments  are  obtained. 

Rewriting 

Nj(t)  =  hni{t)+^fXXj(t)  (4.3) 

Vy(()  =  tey{()+VXy/{f)  (4.4) 

An  approximate  variance  of  the  estimator  of  the  number  of  units  on  node  i. 
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5.  A  Simpler  Poisson  Approximation 

Assume  the  probability  of  allocation  Oi  is  independent  of  iV^t),  Nj(0);  then 
the  number  of  looks  at  node  i  is  a  Poisson  process  with  rate  Aa,-  independent  of 
the  other  nodes.  Further 


ElNiit)]  =  1  +  Oj-f  (5.1) 

V<ir[N,(t)]  =  ai(  (5.3) 

Var[Vi(()]  =  +  o,([iiPi(l  -  Pi)  +  (wf  ]  (5-4) 
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(5.5) 


Co7{Ni{t),Vi{t)]  =  riPiait 


where  we  assume  that  at  time  0  each  node  is  looked  at  once. 

In  this  case  (4.10)  becomes 

Pi  ^  ‘ 

and  the  estimators  fi{t)  are  independent. 

Assume  purposeful  allocation  is  adapted  with  function  hi  as  in  (3.1).  A  simple 

approximation  to  Var  can  be  obtained  by  neglecting  all  covariances  and 

.  » 

assuming  the  number  of  looks  at  node  i,  {Ni(ty,  t  >  0),  is  a  Poisson  process  with 
rate  Aa,-  where  a,-  is  determined  by  (3.6).  In  this  case 

(5.7) 


6.  Numerical  Examples 

Suppose  there  are  3  nodes  with  rj  tmits  on  node  /  with  ri  =  49,  r2  =  25,  and 
r3  =  16.  The  probabilities  of  detecting  a  unit  on  node  /,  pj,  are  p\  =  1/11,  p2  =  0-5, 
andp3  =  10/11. 

The  variance  of  the  estimate  of  the  sum  of  the  units  on  all  the  nodes  under 
purposeful  allocation  with  hi  as  in  (3.1)  was  studied  using  simulation  for  y=  0, 1, 
10.  Each  replication  of  the  simulation  begins  with  one  observation  at  each  node. 
The  times  of  arrival  of  a  Poisson  process  with  rate  1  are  then  simulated.  A  node 
for  observation  at  time  t  is  randomly  chosen  using  probabilities 


ai{t)  =  m 


Oj  ,  Vim- Pi) 
Ni{t)  Niitfpf  _ 
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for  i  =  1,  2,  3  with  K{t)  the  normalizing  constant.  A  binomial  observation  is 

generated  for  the  node  chosen.  The  simulation  has  500  replications. 

Table  1  records  the  sample  mean  and  square  root  of  the  sample  variance  of 

3 

the  svun  of  the  estimates  for  f  =  5, 10, 20, 50  where 

1=1 


Viit) 
m)Pi  ■ 


Also  shovm  are  the  sample  means  of  the  square 
approximate  variances  (1.3) 


Vim-Pi) 
Ni(ifpf  ■ 


root  of  the  sum 


of  estimated 


Table  1 

Estimate  of  Total  Number  of  Units  on  Ail  Nodes 
n  =  49,  rz  =  25,  rs  =  16;  pi  =  1/11,  pz  =  0.5,  ps  =  10/11 


Simulation 

Approximation 

im 

i 

1 

Square  Root  of 
IDifferential 
Equation 
Approx  Variance 

Square  Root 
of  Poisson 
Approx 
Variance 

Time 

1 

Sample 

Mean 

Square  Root 
of  Sample 
Variance 

Sample 

Mean 

5 

89.8 

1 

1 

88.1 

1 

10 

89.3 

15.2 

1 

1^1 

12.5 

10.4 

10 

0 

90.3 

12.4 

11.9 

10.9 

10.9 

1 

89.8 

9.2 

8.2 

9.1 

8.1 

10 

90.1 

10.2 

8.6 

10.2 

8.0 

20 

0 

90.0 

8.9 

8.6 

8.2 

8.2 

1 

89.5 

7.1 

6.1 

6.7 

6.0 

10 

89.5 

8.3 

6.5 

8.0 

6.1 

50 

0 

89.9 

5.8 

5.5 

5.4 

5.4 

1 

89.9 

3.9 

3.9 

4.2 

3.9 

10 

89.5 

6.2 

4.3 

5.7 

4.2 

11 


Table  1  also  displays  results  of  the  approximate  variance  Var 

-  I 

obtained  using  (4.10)  and  (4.12)  and  solving  the  differential  equations  (2.8)  and 
the  second  moment  equations  in  Appendix  A,  (A.27)  -  (A. 32).  The  equations 
were  solved  numerically  using  the  4*^/5*^  order  Runge-Kutta-Fehlberg  method 
as  implemented  in  MATLAB;  (cf.  Math  Works,  1992). 

Table  1  also  displays  the  simple  Poisson  approximation  to  Var  of 

.  i 

(5.7)  with  the  ai  in  (3.6). 

Recall  that  the  estimate  Var  is  unbiased  for  both  approximations.  The 

.  i 

true  value  of  is  90.  The  differential  equation  approximation  is  close  to  the 
square  root  of  the  sample  variance  for  all  values  of  /for  times  10, 20,  and  50.  The 
Poisson  approximation  is  close  to  the  square  root  of  the  sum  of  the  estimated 
approximate  variances;  both  of  these  approximations  are  neglecting  the 
covariances  induced  by  the  purposeful  allocation;  these  covariances  become 
more  pronounced  as  y  becomes  larger.  There  is  no  covariance  for  y=  0,  equally 
likely  allocation.  Note  that  the  Poisson  approximation  is  conservative.  However 
for  y=  1,  the  Poisson  approximation  is  within  about  10%  of  the  differential 
equation  approximation  which  incorporates  the  covariances.  However  the 
difference  is  larger  for  y-  10. 

Note  that  the  square  root  of  the  sample  variance  and  the  differential  equation 
approximation  suggest  that  purposeful  allocation  with  y  =  1  yields  the  smallest 
variance  of  the  estimated  sum  of  the  numbers  of  units  on  all  the  nodes.  A 
rationale  for  this  suggestion  follows. 

Suppose  there  are  a  fixed  number  of  looks  K  that  the  sensor  can  take  of  all 
nodes  and  the  number  of  units  on  each  node  i,  ri  is  known  along  with  the 
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probability  of  detecting  a  unit  on  node  i,  pi.  Let  kj  be  the  number  of  looks  the 
sensor  gives  to  node  j.  If  each  observation  has  a  binomial  distribution 


Var\ 


h(t) 

,i=l 


J  »•=! 


S  hvi 


Lagrange  multipliers  can  be  used  to  show  that  the  (approximate)  ki,  j  =  1, 2, 3 
that  minimize  Var 


im 

Li>l 


are  solution  corresponds  to 


the  a,-  of  (3.6)  with  a,-  =  0  and  y=  1.  Thus,  if  one  is  interested  in  minimizing  the 

estimated  variance  of  the  sum  of  the  number  of  uiuts  on  all  the  nodes,  then  one 
should  look  at  node  i  a  number  of  times  proportional  to  [ri{l-pi)/pif  '^-  If  one 


were  interested  in  minimizing  the  estimated  variance  of  the  estimate  of  the 
number  of  units  on  the  node  with  the  greatest  number  of  units  then  one  would 
allocate  all  looks  to  that  node;  this  corresponds  to  the  purposeful  allocation 
policy  of  y=  «>. 

Tables  2  and  3  present  results  of  the  simulation  experiment  with  ri  =  49, 
r2  =  25,  ra  =  16  and  p\  =  0.7,  p2  =  0.8,  and  pa  =  0.9.  Table  2  presents  the  simulation 
and  approximation  results  for  the  estimate  of  the  sum  of  units  on  all  the  nodes. 
Table  3  presents  the  simulation  and  approximation  results  for  the  number  of 
units  on  the  individual  nodes.  The  differential  equation  results  are  close  to  the 
simulated  values  for  times  f  =  10, 20, 50.  The  Poisson  approximation  also  seems 
to  be  adequate.  The  Poisson  approximation  may  be  doing  better  in  this  case 
because  the  probabilities  of  unit  detection  are  larger.  One  source  of  the 
covariance  between  the  estimators  ^  (f)  is  the  possibility  that  V(i,  t)  may  be  0,  in 
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which  case  node  i  will  not  be  visited  very  frequently  for  y>  0  for  the  purposeful 
allocation  with  function  hi  as  in  (3.1). 


Table  2 

Estimate  of  Total  Number  of  Units  on  All  Nodes 
n  =  49,  rz  =  25,  rs  =  16;  pi  =  0.7,  pz  =  0.8,  pa  =  0.9 


Simulation 

Approximation 

im 

i 

1 

Square  Root  of 
Differential 

Square  Root 
of  Poisson 

Time 

1 

Sample 

Mean 

Square  Root 
of  Sample 
Variance 

Sample 

Mean 

Equation 
Approx  Variance 

Approx 

Variance 

5 

0 

89.9 

5.52 

1 

3.30 

3.30 

1 

89.4 

6.53 

1 

^^9 

3.11 

3.02 

10 

88.9 

8.55 

3.14 

3.06 

2.98 

10 

0 

90.1 

2.75 

2.84 

2.59 

2.59 

1 

89.9 

2.37 

2.45 

2.46 

2.35 

10 

89.7 

2.46 

2.53 

2.51 

2.35 

20 

0 

90.1 

2.04 

2.07 

1.95 

1.95 

1 

90.0 

1.72 

1.80 

1.90 

1.76 

10 

89.9 

1.87 

1.88 

2.06 

1.80 

50 

0 

90.0 

1.30 

1 

1.28 

1.28 

1 

90.0 

1.12 

1 

mm 

1.29 

1.16 

10 

89.9 

1.26 

L 

1.24 

1.35 

1.21 
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50 

20 

10 

CJl 

Time 

■ 

o 

10 

■ 

o 

10 

B 

o 

10 

t-l 

o 

49.0 

49.0 

49.1 

49.0 

49.0 

49.1 

48.9 

49.0 

49.1 

48.5 

48.6 

49.0 

Mean 

' 

\ 

O  O 

0.88 

(0.85) 

i-a 

o  b 

VO  00 

H-i 

3  ON 

b  b 
o  ro 

(99T) 

89T 

In  b 
3  ^ 

"C  i-k 

^  2 
^  h-a 

ro  ro 
ro  b 
o  w 

4.85 

(2.00) 

3.95 

(2.27) 

3.89 

(2.81) 

Square  Root 
Sample 
Variance 
(Di£f  Eqtn 
Approx) 

o  o 

3^ 

0.86 

(0.86)__ 

H-k 

b 

o  ro 

3  ^ 

1—1 

b  b 
w  ro 

♦— »  l-i 

On  (-n 

"C  1—1 
b  b 

OJ  t-i 

'C  t— 1 
bo  bo 

H-l  H-l 

1o  ro 

b5  ^ 

O  VO 

'w' 

2.12 

(2.17) 

2.37 

(2.38) 

3.11 

(2.81) 

Sample 

Mean 

(Poisson 

Approx) 

~.(sj 

24.9 

25.0 

25.0 

24.9 

24.9 

25.0 

24.9 

24.9 

25.0 

24.6 

bo 

24.9 

Mean 

*::t> 

O  O 

O  O 
ON  On 

3  p 
b  b 
3  w 

Cl  1—1 

b  b 

w  ^ 

P  p 

b  b 

On  vq 

O  O 

b  b 

O  00 

b 
^  o 

C*  »-i 
ro  b 

ON 

C»  ^ 
ro  b 

O  ON 

2.94 

(1.88) 

ro 
b  b 
3  00 

2.20 

(1.53) 

Square  Root 
Sample 
Variance 
(Diff  Eqtn 
Approx) 

o  o 

o  ro 

o  o 

ON  ON 

o  p 

In  b 
\o  *-» 

'C  H-* 

b  U 
CJl  o 

0.98 

(0.95) 

o  o 
b  b 

O  Ol 

(—1 
b  4i>^ 
3  00 

H-l  1-1 

ro  b 
b  ro 

C  1—1 

ro  b 
o  o 

(OZT) 

68T 

1.69 

(1.59) 

(€fiT) 

89T 

Sample 

Mean 

(Poisson 

Approx) 

_ 

Q> 

~-N> 

16.0 

16.0 

16.0 

15.9 

16.0 

16.0 

15.9 

16.0 

16.0 

15.8 

16.0 

16.0 

1 

g 

0.71 

(0.66) 

0.45 

(0-45) _ 

0.33 

(0.32) 

1.06 

(1.00) 

0.71 

(0.68) 

0.52 

(0.48) 

1.33 

(1.29) 

0.96 

(0.89) 

0.71 

(0.64) 

1.97 

(1.33) 

1.54 

(1.07) 

1.17 

(0.82) 

Square  Root 
Sample 
Variance 
(Diff  Eqtn 
Approx) 

0.65 

(0.61) 

0.46 

(0.44) 

0.32 

(0.32) 

0.95 

(0.84) 

0.71 

(0.65) 

0.50 

(0.48) 

1.26 

(1.00) 

0.93 

(0.82) 

0.68 

(0.64) 

1.32 

(1.13) 

(66*0) 

€IT 

0.87 

(0.82) 

Sample 

Mean 

(Poisson 

Approx) 
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Table  3 

Simulation:  500  replications 

_ Estimate  of  Number  of  Units  on  Node _ 

Nodel  I  Node!  _ Node  3 
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APPENDIX  A 


In  the  Appendix  we  present  details  of  the  normal  approximation.  We  follow 
an  analytical  approach  used  in  a  different  context  in  Gaver  and  Jacobs  (1995),  and 
in  Gaver,  Morrison,  and  Silveira  (1993).  Let  the  moment-generating  function 
(assumed  to  exist,  otherwise  use  the  characteristic  function)  be 


Vf(fte;()  =  E[exp{6N(()+?V(f)}] 


■ 

\l  1  1 

E  exp 

iww+swo 

[/=1  ;=1  J. 

Condition  on  (N,(f),  W(0),  i  e  {1, 2, . . .,  1}  to  obtain 


1 

exp^ 

X«;Ny{f  +  h)  +  {,Vy((  +  fl) 

[pl  J 

=  (l-A/i)exp{fflV(f)  +  ^(f)}  (A.2) 

+AAXo.(NW.V(<))exp{eN(()+?V(f)}[e*‘t(6)] 

i 

where 

=  (A.3) 

with  Z(i)  an  observation  of  the  number  of  units  on  node  i. 

Let  ft  -4  0  to  obtain 


Scaling 

Let 


(A.5) 
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and  let  A»  1. 


Let 

p(««;()  =  E[exp{«X(()+W)}]; 

then 

\if(el  VI,  ;  f )  =  ((>{e,  t)exp{VI[ftw(f)  +  ^(t)]}  • 

Thus,  we  have  the  following  equation  from  (A.8)  and  (A.4) 


=  ^<P{B.i:t)exp[4l[em{l)  +  {!)(()]} 

+^(p{0,  $;t)exp{VI[ft«(0  +  $y(f)]}[ftM'(t)  + 


(A.6) 

(A.7) 

(A.8) 


xex 

Dividing  both  sides  by  exp|VI[ft«(f)  +  5^(0]}  obtain 

^(p{0,  ) + VI  (p{0,  t)[0m'{t)  +  ^'(0] 

1=1  [[  j  (A.IO) 

+^^a,(m(f),n(f))^Yy(f)+0^^j  -exp{fiX’(f)+^(f)}  e®‘^4(^,7VI) 

j  ^ 
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Let 


ai{m{t),v{t))  = 


(A.ll) 


where  hfx,  y)  is  a  sufficiently  smooth  function  with  first  order  partial  derivatives 
and  hj^^x, y j  =  XFhj [x, y) . 


Hi{x;t)  = 


■£^hi(mi{t),Vi{t)) 


(A.12) 


(A.13) 


(A.14) 


(A.15) 


Note  that 


a,(Am,(f) + 4XXi{t),Xvi{t) + ^[XYi{t)) 

hi{Tni{t),Vi{t)) 

;=i 


(A.16) 


h  Hi{x;t)Xi{t)+Hi(y;t)Yi{t)+J^Hij,{x;t)Xk^^^^ 
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Since 


20 


37  ^o(®/  + ^"(0] 

at 

=  -X<po{e,^;t) 

+ ~^i^i  ■•■  ) + 2x(^‘^  ■•■  ) 

Z^K(0/^;(0) 


(A.20) 

Equating  terms  of  order  A^/2,  the  terms  of  order  A  cancel.  The  terms  of  order  ^[X 


result  in  the  equation 


(po{O,^}t)l0nt'{t)+^Xt)] 

f=i  X^/K(0/i'y(0) 


*Hi{x-.t)^(i>om,t]-<-iii{r,t)-^(i>o{e,i;t) 

oui  dt^i 

+i  %(x;()^«>o(««;0+Htt(y;')^«)(»-6')  •+o(x} 

* 


(A.21) 
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Equations  for  the  joint  moments  of  [iXft),Yj(t)}  can  be  obtained  by 
differentiating  (A.26)  with  respect  to  [Oi)  and  and  evaluated  at  6=^=0.  The 
resulting  equations  are 
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(A.27) 


h(mi{t),vg{t)) 


+  2H^(a:;0£[x^(0^]  +  2H^(y;0E[^^(O>'KO] 


7 

+2Z««(^:')E[X<(()Xt(()]+2SHtt(y;f)E[X/(()>i(»)] 

it  It 


For  i*a 


|e[X<(()X.(()] 

=  H„(A:;()E[X.(f)XX()]+H.(y;<)£[n(<)X<(()] 
+H<(j:;()E[X,(()X<(()]+H,(y;()E[X„(()Y<(()] 

+lHfl)c(^;')E[Xt(t)X<(()]+lH,i(y;()E[Xi(()Vt(()] 
k  k 

+lHA(>;;0£[X«(()Xn(()]  +  SHtt(y;()E[X,(f)yt(()]. 

it  It 


E(m, ((>,«))  ) 

XE(my(f),i>y(()) 


+2i,(E)[H<(i;()E[X,(()r<(()]  +  H,(y;/)E[y7(()]] 

+2i>,(<)XH<y(Ar;()E[x,(f)y,(()]+H</(y;()£[vi(()V,(f)] 

7 


(A.28) 


(A.29) 
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Tot  j^k 


=  {H,(;(;()£[Xy(f)yt(f)]  +  Hy(y;()£[y,(/)yt(f)]}i>iO) 
+{Ht(jr;()£[Xt(()y,(<)]  +  Ht(y;l)E[y;(t)y|iW]}f>lW 
+*>i(/)S«>((>^;‘)E[x<(f)yt(0]+«,7(y;')E['^t(')n(‘)] 

t 

+i^WlHm(:t;()E[x<(0'',(<)]+Ht,(y;f)E[y,(()y;(()] 

t 

j 

+Ht(>:;(){E[yt(()Xt(()]  +  £[x|(()]f>,(t)} 

+Ht(y;  ()[£[n^(<)]  +  E[Xt(()yi(/)]l>i(*:)] 

+X{%(*;<)E[Xy{()yt(<)]  +  Hiy(y;0£[yt(()y, ■(()]} 

/ 

+l{Hi,{Ar;f)£[Xy(f)Xt(()]  +  %(y;()E[Xt(f)y,(()]]i^(E) 


(A.30) 


(A.31) 


24 


For 


^£[X,(()V|t(()] 

=  H<(l;/)E[X<(l)Yt(()]  +  H,(y;()E[y,(f)yt(()] 
+{H*(i;()E[X<(()Xn(l)]  +  Ht(y;OE[X,(()yt(()]}f>l{*:) 

+XH<y(j(;()E[X;(t)y»(()]  +  H</(y;f)Eln(<)V,(<)] 

;■ 

+EK(*;0E[X/(0X<(0]+%(y;')E['^;(')X/{0]]il(*:) 


(A.32) 
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